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ABSTRACT 

We present a Keck II LRIS spectroscopic follow-up study of the possible optical counterparts to a 
flux-limited sample of galaxies selected from an 850-/mi survey of massive lensing clusters using the 
SCUBA bolometer array on the JCMT. These sources represent a population of luminous dusty galaxies 
responsible for the bulk of the 850-^m background detected by COBE and thus for a substantial fraction 
of the total far-infrared emission in the Universe. We present reliable redshifts for 20 galaxies and redshift 
limits for a further four galaxies selected from the error-boxes of 14 submillimeter (submm) sources. Two 
other submm detections in the sample have no obvious optical counterparts, and the final submm source 
was only identified from imaging data after the completion of our spectroscopic observations. The 
optical identifications for 4 of the submm sources have been confirmed through either their detection 
in CO at mm-wavelengths (two pairs of galaxies at z = 2.55 and z — 2.80) or from the characteristics 
of their spectral energy distributions (two of the central cD galaxies in the lensing clusters). Plausible 
arguments based on the optical spectral properties (starburst or AGN signatures) of the counterparts 
allow us to identify a further two likely counterparts at z = 1.06 and 1.16. For the remaining 8 cases, 
it is not always clear which, if any, of the optical sources identified are the true counterparts. Possible 
counterparts for these have redshifts ranging from z = 0.18 to z = 2.11. The application of a range of 
techniques, including near- and mid-infrared imaging and radio mapping, will assist in the identification 
of the true sources of the submm emission, while CO line mapping with current mm-interferometers and 
hard X-ray observations should aid in the determination of the nature of their emission. Working with 
the current identifications, we suggest that the majority of the extragalactic background light in the 
submm is emitted by sources at z < 3 and hence that the peak activity in highly-obscured sources (both 
AGN and starbursts) lies at relatively modest redshifts. We find that a lower limit of 20 per cent of the 
submm sources in our sample show some sign of AGN activity; however, we caution that this does not 
necessarily translate into a 20 per cent AGN contribution to the measured submm emission from these 
sources. 

Subject headings: cosmology: observations — galaxies: distances and redshifts — galaxies: evolution — 
galaxies: formation — galaxies: active — galaxies: starburst 



1. INTRODUCTION 

The cumulative emission from all objects lying beyond 
the Galaxy, the extragalactic background light (EBL), pro- 
vides important constraints on the integrated star forma- 
tion history of the Universe. The recent measurement of 
the EBL at far-infrared (FIR) and submm wavelengths us- 
ing data from the FIRAS and DIRBE experiments on the 
COBE satellite (Puget et al. 1996; Guiderdoni et al. 1997; 
Schlegel et al. 1998; Fixsen et al. 1998; Hauser et al. 1998) 
indicates that the total emission from star formation and 
AGN activity that is absorbed by dust and reradiated into 
the FIR/submm is comparable to the unobscured emission 
seen in the optical. This suggests that obscured star for- 
mation may be responsible for a large fraction of the stars 
and metals seen in the local Universe, a conclusion which 
would have profound consequences for models of galaxy 
formation and evolution. 



Independent support for the claim that optical estimates 
of the star formation density in the distant Universe may 
be missing a substantial component that is obscured by 
dust comes from deep submm surveys with the new cam- 
era SCUBA (Submillimeter Common User Bolometer Ar- 
ray; Holland et al. 1999) on the 15-m James Clerk Maxwell 
Telescope 7 (JCMT) on Mauna Kea. SCUBA has, for the 
first time, enabled deep, unbiased surveys to be made of 
the submm sky (Smail, Ivison & Blain 1997; Barger et al. 
1998, 1999; Hughes et al. 1998; Eales et al. 1999). These 
surveys have uncovered numerous sources with proper- 
ties similar to those expected for distant ultraluminous 
infrared galaxies (Barger et al. 1998; Smail et al. 1998). 
If the majority of the submm emission in these systems 
comes from dust-obscured star formation, then their in- 
ferred star formation rates are of the order of several hun- 
dred solar masses per year. 
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Fig. 1. — 15" X 15" images of the 16 submm sources from the Smail et al. (1998) cluster survey that have corresponding optical imaging. 
The orientations of the images are arbitrary. The objects are ordered from the upper-left by decreasing apparent 850-fim flux. Error circles of 
3" radius are centered on the submm source positions. The possible optical counterparts in each field are labeled and can be cross-referenced 
with the data in Table 2. 



The deepest submm counts (Blain et al. 1999b), fluctu- 
ation analyses (e.g. Hughes et al. 1998), and analyses of 
deep and wide-area surveys (e.g. Barger et al. 1999) indi- 
cate that the bulk of the background emission at 850-/im 
detected by COBE is resolved into discrete sources at a 
flux limit of < 1 mjy. Thus, we are in a position to un- 
dertake detailed studies of the population responsible for 
the majority of the emission in the FIR background. In 
particular, we can measure the redshift distribution of the 
submm population and use this to trace the extent and 
evolution of obscured star formation in the distant Uni- 
verse. The detailed study of the resolved component of 
the background also provides a clear method for deter- 
mining what fraction of the submm emission originates 
from AGN activity rather than star formation (Blain et 
al. 1999a; Guiderdoni et al. 1998). 

Finally, it is conceivable that some fraction of the popu- 
lation responsible for the FIR background lies at high red- 
shift, z ~ 5-10. Submm selection is a powerful technique 
for locating such distant galaxies. Indeed, the steep ther- 
mal dust spectrum, which peaks in the FIR at a rest-frame 
wavelength of about 100-/mi, is redshifted into the submm 
for z > 1. The resulting strong negative K-correction for 
sources out to z ~ 10 is sufficient to offset cosmological 
dimming for qo = 0.5. Even for low values of qo the 850- 
/zm flux density is only expected to decrease by a factor 
of a few over this redshift range (Blain & Longair 1993; 
Hughes, Dunlop & Rawlings 1997). 

In this paper we present spectroscopic redshift infor- 
mation for a sample of submm-selected galaxies from the 
SCUBA cluster lens survey of Smail et al. (1997, 1998). 
One advantage of a lensed survey is that the problem of 
source confusion (which, due to the coarse resolution of 
submm telescopes, can contribute noise in faint images; 
Blain, Ivison, & Smail 1998) is reduced since both the flux 
densities and the mean separations on the sky of the back- 



ground sources are increased (Blain et al. 1999b). The 
primary advantage, however, is that the clusters magnify 
any background sources (here the median amplification of 
the submm fluxes is ~ 2.5), thereby providing otherwise 
unachievable sensitivity in the submm and easing spectro- 
scopic follow-up in the optical. The full 850-/im survey 
detected 17 sources above 3c significance (la < 2 mjy in 
the image plane) over a total surveyed area of 36arcmin 2 . 
Crude redshift limits for these sources were derived from 
broad-band imaging by Smail et al. (1998), who inferred 
that most of the galaxies were at z < 5. 

Section 2 describes the submm sample and the spectro- 
scopic follow-up. Section 3 presents the redshift identifi- 
cation for each submm source. Finally, section 4 discuss 
the results and gives our conclusions. 

2. SAMPLE AND OBSERVATIONS 

Smail et al. (1998) sought the optical counterparts for 
their 17 submm detections within a conservative <J 6 arcsec 
radius (the combined 2a positional uncertainty of their 
faintest sources) from the nominal submm position. All 
possible optical counterparts were identified in deep HST 
and ground-based imaging data to / ~ 26.0 and 23.5, 
respectively, for the 16 submm sources covered by ex- 
isting optical imaging data. Using a deep Keck II LRIS 
/-band image, we recently identified a probable opti- 
cal counterpart (I = 24.3) to the seventeenth source, 
SMM J0443+0210 (N5; Fig. |l|). Unfortunately, this iden- 
tification came too late for the source to be included in 
the spectroscopic campaign described here. Two of the 
17 sources have no obvious counterparts to the magnitude 
limits of their respective images. 

Since the lensing clusters lie outside the Galactic plane, 
the 850-/im sources are unlikely to arise from local regions, 
and indeed we find no evidence for dark patches caused 
by dust foreground globules in any of the optical images. 
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Table 1 

Log of the Keck II LRIS spectroscopic 

OBSERVATIONS 



UT Date 


Seeing 
(") 


Targets 


Exposure 
(ks) 


1998 July 18 


0.6 


J1/J2, ,73 


3.6 






K2, K3 


3.6 






PI, P2 


3.6 






M4, M6, M10 


1.8 






Ml, M2 


2.4 


1998 Aug 21 


0.6 


PI, P4 


3.6 






L3, L4 


3.6 


1998 Aug 22 


0.6 


PI, P2 


3.6 






Ml, M3, Mil 


2.4 


1998 Sept 16 


0.6 


PI, P2 


3.1 






Ml, M3 


3.6 


1998 Sept 17 


0.6 


K2, cD 


3.6 






Ml, M3 


3.6 






Nl. N2 


3.6 


1998 Oct 22 


0.6 


H4 


3.0 



Furthermore, the weak or non-detection of these sources 
at 450-/xm (Smail et al. 1997) rules out their origins being 
local. 

As noted by Smail et al. (1998), since the optical coun- 
terparts to the submm sources are likely to be faint and 
the cluster fields are crowded, there is some likelihood that 
unrelated galaxies will fall within the submm error boxes. 
The probability of a galaxy with a given observed apparent 
magnitude m falling at random within a circle of radius r 
centered on the submm source is proportional to r 2 times 
the number density per surface area on the sky of galaxies 
with magnitude m. Smail et al. estimated on an object- 
by-object basis (from the observed galaxy number counts 
for each optical image) the probability that a galaxy with 
the observed apparent magnitude or brighter would fall at 
random within a circle defined by the optical and submm 
positions. Although this procedure provides some measure 
of the reliability of the identifications, it does not allow the 
definitive determination of the true counterparts to indi- 
vidual submm sources. Thus, we have simply used these 
probability estimates to determine the priority with which 
to target sources, but in most cases we have targeted all 
visible counterparts within the error circles. 

The optical spectroscopy of the candidate submm 
sources was undertaken with the Low-Resolution Imag- 
ing Spectrometer (LRIS; Oke et al. 1995) on the Keck II 
10-m telescope using a wide 1.5 arcsec long-slit during sev- 
eral runs from 1998 July to 1998 October. With the wide 
slit and the 300 lines mm -1 grating blazed at 5000 A, the 
resolution was 14 A. The wavelength coverage was varied 
slightly, depending on the target and whether there were 
additional emission lines to search for. Two objects were 
observed per slit, which defined the position angle. The 
observations were typically lhr per slit, broken into three 
sets of exposures. Some of the objects were re-observed 
several times to improve the signal-to-noise. Table 1 gives 
a log of the Keck spectroscopic observations, including 
date, seeing, targets, and exposure times. The objects 
were stepped along the slit by 10 arcsec in each direction, 
and the sky backgrounds were removed using the median 
of the images to avoid the difficult and time-consuming 
problems of flat-fielding LRIS data. Details of the spec- 



troscopic reduction procedures can be found in Cowie et 
al. (1996). 

Figure^ shows the fields for the 17 submm sources from 
the SCUBA cluster lens survey, including the newly identi- 
fied source SMM J0443+0210. The individual images are 
centered on the submm positions with nominal 3 arcsec 
(1<t) radius error circles marked. Each square box covers 
an area 15" x 15"; the orientations of the images are ar- 
bitrary. The possible optical counterparts for each of the 
submm detections are labeled by name and can be cross- 
referenced with Table 0. In both Fig. |l] and Table || the 
submm sources are ordered by apparent 850-/im flux for 
the 4(T and 3a samples. The columns in Table |^ include 
the submm source name, apparent 850- /im flux, candidate 
optical counterparts (most likely listed first), / magnitudes 
for the counterparts, redshifts for the counterparts, name 
and redshift of the cluster field, and lensing amplifications 
for the counterparts at the given redshifts. 

3. REDSHIFT IDENTIFICATION 

Our ability to conclude whether a particular galaxy is 
likely to be the submm source depends strongly on whether 
the optical spectrum of the galaxy shows any remarkable 
features, such as particularly strong [On] A3727 or Ha 
emission lines that indicate a starburst, or high excita- 
tion or broad lines that show the presence of an AGN. 
The relative paucity of AGNs in the general field popu- 
lation (<; 1%) suggests that if an AGN is identified, then 
the submm emission is most likely associated with that 
source. We note, however, that the visibility of remarkable 
features, especially AGN lines, will be strongly affected 
by dust in the galaxy - the same component which our 
submm selection should guarantee is present. The iden- 
tification process is therefore somewhat problematic, and 
while we are capable of robustly identifying some submm 
sources with optical counterparts, specifically those with 
the most striking spectral features (as has been confirmed 
through the CO detection of two of our candidates), this 
approach leaves us with a number of ambiguous cases. As 
a secondary criterion, we consider the morphologies of the 
objects and the presence of any merger activity, which can 
be an indicator of luminous FIR systems in the local Uni- 
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Fig. 2. — Keck spectra, where available, of the most likely optical counterparts to the submm detections: (a) M2 (z = 1.23), (b) J1/J2 
(z = 2.55), (c) L3 (z = 1.06), (d) P4 (2 = 1.16), (e) Nl (z = 0.18), (f) K2 (z = 1.60?), (g) M6/M10 (z = 0.21), and (h) P2 (z = 2.11?). 
For (b) J1/J2 and (g) M6/M10 the spectra of both members of the galaxy pairs are shown. The spectrum for J3 is not shown because the 
redshift is inconclusive. 



verse. 

We now discuss the nature and reliability of our redshift 
measurements in order of decreasing submm flux density, 
as presented in Fig. [l] and Table |[ In Fig. || we show the 
Keck spectra for the most likely optical counterparts to 
the submm sources, as discussed below. 

SMMJ02399-0136: A redshift of z = 2.80 was mea- 
sured for the pair of optical sources, L1/L2, by Ivison et 
al. (1998). LI is a compact, dust-obscured AGN, and L2 



is a companion structure whose extended emission may 
result from a vigorous burst of star formation triggered 
by an interaction with LI. The detection of CO emission 
in the mm-waveband coincident in redshift and position 
with the optical counterpart (Frayer et al. 1998) leaves no 
doubt that the L1/L2 pair is the correct identification. 
SMM J09429+4658: The redshift for the nearby bright 
galaxy HI is z = 0.33 (Dressier et al. 1999). It shows no 
[On] A3727 or H a emission but contains an obvious dust 
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Fig. 3. — The 2D spectral image for M3 (z = 0.94) with tick marks indicating the positions of [On] A3727, H/3, and [Oin] A5007. Although 
H/3 cannot be reliably identified, the emission line identified with [Oil] is unambiguous and that identified with [Om] is also very clear. 



lane and thus might plausibly be the submm source. H4 
shows Ha and NaD in absorption and is a blue star. An- 
other plausible counterpart — a faint red source close to 
the nominal submm position — has recently been uncov- 
ered in deep near-infrared imaging of the SCUBA fields 
(Smail et al. 1999). Given the dust characteristics of HI 
and the very low probability of seeing such a bright fore- 
ground galaxy close to the submm position, we consider 
HI to be the most likely counterpart; however, we caution 
that there is some uncertainty about this identification. 
SMM J14009+0252: We targeted J3, the nearest galaxy 
to the nominal submm position, and found the spectrum 
to be featureless and flat, so no redshift identification could 
be made; however, from its blue SED and lack of obvious 
emission lines, the object does not appear to lie at z -C 1.5 
or z 3> 2.5. No spectroscopic observations were made of 
J4. 

SMM J14011+0252: We identify this submm source 
with the high redshift starburst pair J1/J2 at z = 2.55, 
based on strong Lya emission and both CivA1548 and 
Si IV A1396 absorption features (Fig. |^b) . This pair of clas- 
sic Lyman break galaxies was subsequently detected in CO 
at close to the optical redshift (Frayer et al. 1999). As with 
SMM J02399— 0136, the spatial and redshift coincidence of 
the CO emission confirms our identification. 
SMM J02399 0134: The galaxy L3 at a redshift of 
z = 1.06 is almost certainly the submm source. The 
galaxy has an unusual ring morphology (Fig. [l]), and the 
nucleus hosts a Seyfert 1.5 AGN, characterized by strong 
[O n] A3727, [Ne v] A3426, and [Ne III] A3869 emission, in 
addition to strong Mgn A2800 absorption (Fig. ||c). L4 is 
at z = 0.42 and hence is background to the cluster. L5 is 
a passive cluster elliptical at z = 0.37 (ID #32 in Mellier 
et al. 1988). 

SMM J22471 0206: P4 is a z = 1.16 emission-line 
galaxy with [O n] A3727 and Mg n A2800 in absorption over 
weak broad Mgn emission (Fig. ||d). The highly peculiar 
morphology of this source, plus its weak AGN character- 
istics, suggest that it is the correct optical counterpart. 
Other candidates include PI, which has a flat featureless 
spectrum, probably putting it at z ~ 2, and the more dis- 
tant (3.3") and less probable P3, which was not observed. 
SMM J21536+1741: This source lies close to the cen- 
tral cD galaxy of the cluster A 2390 (z = 0.23). Since the 
cD galaxy is a strong radio source, there is little doubt 
that this is the correct identification. The submm emis- 
sion from this source and its relation to the overall spectral 
energy distribution of the central galaxy are discussed in 



more detail in Edge et al. (1999). 

SMM J02400— 0134: There is no obvious optical coun- 
terpart for this source down to the deep limit of the HST 
image. 

SMM J04431+0210: The bright spiral galaxy Nl shows 
strong H a and [O n] A3727 emission (Fig. 0e) and is a clus- 
ter member at z = 0.18 (see also Gioia et al. 1998). The 
probability of seeing such a bright cluster galaxy only 2.3" 
from the nominal submm position is very low. However, 
this is another field where the deep near-infrared imaging 
by Smail et al. (1999) has uncovered a faint red source 
within the submm error-box. Thus, we caution that there 
is still a question mark over the identification of Nl with 
the submm emission. 

SMM J21536+1742: K2 has a relatively uncertain red- 
shift identification of z = 1.60 based on weak HenA1640 
emission and on AlmA1853 - 1862 and MgnA2800 ab- 
sorption features (Fig. &) . There are no outstanding spec- 
tral characteristics marking this object as the true opti- 
cal counterpart. K3 is a clear absorption-line galaxy at 
z = 1.02, which is also detected in the deep ISO CAM 
observations of this field (J. -P. Kneib, unpublished). How- 
ever, its large separation (6.3") from the submm position 
makes it a less likely counterpart. 

SMM J00265+1710: M6/M10 are a pair of emission- 
line galaxies at z = 0.21 (Fig. ||g), foreground to the clus- 
ter, which appear to be interacting. We were unable to 
obtain a redshift for the other candidate counterpart, M7, 
because of its close proximity to M6. The apparent in- 
teraction of the M6/M10 pair may be sufficient justifica- 
tion to suggest that this is the optical counterpart to the 
submm source. However, the spectrum of neither galaxy 
is particularly unusual, and we note that the galaxies are 
not associated with strong 15-^m emission in the deep ISO 
CAM image of this cluster (J. -P. Kneib, unpublished). If 
the SCUBA source is associated with these galaxies then 
the absence of 15-/im emission is unusual given their low 
redshift. For this reason we caution that, while attractive, 
the identification of the submm emission with M6/M10 
may be incorrect. 

SMM J22472 0206: P2 is the only likely visible optical 
counterpart to the submm detection. Our best estimate 
is that it is a star- forming galaxy at z — 2.11, based on 
what appears to be Lya emission and the Lyman break 
(Fig. ^h). However, there is no absorption- line confirma- 
tion, and the ultraviolet end of the spectrum is obtained 
from an integration of only 1 hr. This redshift identifica- 
tion is therefore uncertain. 
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SMM J00266+1710: M3 is an emission-line galaxy at 
z = 0.94 with [O n] A3727 and [O III] A5007 emission lines. 
There is no sign of Mg n A2800 absorption or emission in 
this galaxy, but we believe that the redshift identification 
is reliable. In Fig. || we show the 2D spectral image for 
this galaxy with tick marks at the positions of [O n] A3727, 
H/3, and [O in] A5007. M3 is the only visible optical coun- 
terpart in the vicinity of the submm position. 
SMM J14014+0252: The core of the cluster targeted 
in this field, Abelll835, contains a massive cooling flow 
centered on the z = 0.25 cD galaxy close to the submm 
source. The strong line-emission and other spectral fea- 
tures of massive star formation detected in this galaxy 
(Allen 1995) suggest that it is the source of the submm 
emission. Edge et al. (1999) identify a cool dust contribu- 
tion in the observed submm spectral energy distribution 
of this galaxy and suggest that it arises from dust heated 
by star formation. 

SMM J00267+1709: There is no visible optical coun- 
terpart to this source on the deep HST F814W exposure, 
placing a limit of / > 26 on the apparent magnitude of 
any counterpart. 

SMM J04433+0210: This is the faintest submm source 
to make it into the catalog. Moreover, this source lies off 
the original HST WFPC2 exposure of this cluster used by 
Smail et al. (1998) to identify candidate counterparts. A 
recent deep Keck /-band image indicates that an / = 24.2 
object lies within 1.2" of the submm position. There are no 
brighter sources within 6" of this position. The probabil- 
ity that this faint object is the counterpart of the submm 
source is P = 0.04 (calculated in the same manner as 
in Smail et al. 1998). This suggests that the faint opti- 
cal source may be associated with the submm emission. 
However, as this candidate was only recently acquired and 
is extremely faint, no spectroscopic observations have yet 
been undertaken. 

Since many of our most probable counterparts do not 
show any unusual spectral features, it is not always cer- 
tain that we have identified the true optical counterpart to 
each submm source. Two of the submm detections in the 
sample have no visible optical counterparts in very deep 
imaging, and it is possible that the true counterparts to 
some other sources in the sample are similarly optically 
faint. In principal such sources could either be at very 
high redshift or be so highly obscured that they are emit- 
ting their energy almost entirely in the submm (Dey et al. 
1999). 

4. DISCUSSION AND CONCLUSIONS 

We have carried out a spectroscopic study of the candi- 
date optical counterparts for 14 of the 17 submm sources 
(excluding the two blank fields and the unobserved coun- 
terpart to SMMJ04433+0210) from a gravitationally- 
lensed submm-selected survey to determine redshifts and 
crucial spectral information for a representative sample of 
the submm source population. This information is essen- 
tial for understanding the evolution of obscured star for- 
mation in the Universe and the contribution of AGN to 
the FIR background. 

Our survey has produced firm identifications for 
four of the submm sources: two central cD galax- 
ies in the lensing clusters (SMM J21536+1741 and 
SMM J14014+252), an interacting pair of galaxies at z = 



2.80 (SMMJ02399-0136 L1/L2), one of which hosts a 
Seyfert-2 nucleus (LI), and a further pair of galaxies 
at z = 2.55 (SMMJ14011+0252 J1/J2) that show star- 
burst features. It has also produced reliable identifica- 
tions for a further two weak AGN sources at z = 1.06 
(SMMJ02399-0134L3) and z = 1.16 (SMM J22471-0206 
P4). In the following we refer to the four non-cluster 
sources listed above as our reliable sample. 
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Fig. 4. — The top panel shows the redshift distribution for the most 
likely optical counterparts to the submm detections in our spectroscopic 
survey. The cross-hatched histogram represents the reliable sample of 
four non-cluster counterparts (L1/L2, J1/J2, L3, and P4) and the open 
histogram represents cluster contamination (Nl and the two cDs). The 
redshift for J3 is uncertain and only expected to be in the range z — 1.5— 
2.5; thus, it has been plotted at z — 2. The three models discussed 
in the text, including Model E from Guiderdoni et al. (1998) (dotted), 
the original Gaussian model from Blain et al. (1999a) (dashed), and the 
revised Gaussian model (solid curve; see text) are superimposed on the 
redshift distribution in panel one. The lower two panels give the re- 
sults published from the SCUBA survey of the Hubble Deep Field by 
Hughes et al. (1998) and from the on-going survey of the CFRS fields 
by Lilly et al. (1999). The shaded histograms in the bottom two panels 
represent the spectroscopic identifications and the dotted histograms the 
photometrically-derived limits. The number of blank fields in our survey 
and in the CFRS survey are indicated by floating boxes at nominally 
high redshift. 



The remaining eight sources have candidate counter- 
parts with redshifts ranging from z = 0.18-2.11, but they 
lack remarkable spectral characteristics to clearly identify 
them as the true counterparts. The lower redshift (z <C 1) 
systems include a pair of interacting galaxies at z — 0.21 
(SMM J00265+1710 M6/M10) and two bright spiral galax- 
ies at z = 0.18 (SMM J04431+0210 Nl; cluster member) 
and z = 0.33 (SMM J09429+4658 HI). 

In the top panel of Fig. |] we show our differential red- 
shift distribution for the most likely optical counterparts 
to our sample of lensed submm detections (excluding the 
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Table 2 

Redshift catalog for the candidate optical counterparts to the submm sample 



Submm 


850 (im 


Candidate Optical 


/ 


Counterpart 


Cluster 


Lensing 


Source 


Flux (mjy) 


Counterparts 


Mag 


Rcdshifts 


Name /Redshift 


Amplification 


4cr Detections 














SMM J02399— 0136 


25.4 


L1/L2 


20.41 


2.80/2.80" 


A 370/0.37 


2.4 


SMM J00266+1708 


18.6 


M2 


22.38 


1.23 


CI 0024+16/0.39 


1.6 






Mil 


22.56 


1.06 


1.5 






Ml 


21.99 


0.39 




1 






M8 


21.99 


0.44 




1.05 


SMM, J 09429+4658 


17.2 


HI 


19.24 


0.33 b 


CI 0939+47/0.40 


1 






H4 


22.15 


star 






SMM J14009+0252 


14.5 


J3 


23.55 


1.5-2.5 


A 1835/0.25 


< 1.7 






J4 


21.62 






SMM J14011+0252 


12.3 


J1/J2 


20.32 


2.55/2.55 


A 1835/0.25 


2.7 


SMM J02399 — 0134 


11.0 


L3 


20.52 


1.06 


A 370/0.37 


2.5 






L4 


21.50 


0.42 


1.1 






L5 


19.00 


0.37 c 




1 


SMM J22471 — 0206 


9.2 


P4 


21.72 


1.16 


CI 2244-02/0.33 


1.9 






PI 


22.92 


~ 2? 


2.3 






P3 


23.17 








SMM J21536+1741 


9.1 


cD 


15.94 


0.23 


A 2390/0.23 


1 


SMM J02400 — 0134 


7.6 




> 26 


z > 4? 


A 370/0.37 


> 1.9 


SMM J04431+0210 


7.2 


Nl 


18.42 


0.18 d 


MS 0440+02/0.19 


1 


3a Detections 














SMM J21536+1742 


6.7 


K2 


24.69 


1.60? 


A 2390/0.23 


1.9 






K3 


21.36 


1.02 


1.7 


SMM .100265+1710 


6.1 


M6/M10 


20.53 


0.21/0.21 = 


CI 0024+16/0.39 


1 






M7 


21.19 








SMM J22472 — 0206 


6.1 


P2 


24.05 


2.11? 


CI 2244-02/0.33 


2.2 


SMM J00266+1710 


5.9 


M3 


23.08 


0.94 


CI 0024+16/0.39 


3.6 


SMM J14010+0252 


5.4 


cD 


15.50 


0.25 


A 1835/0.25 


1 


SMM J00267+1709 


5.0 




> 25 


z > 4? 


CI 0024+16/0.39 


> 2.2 


SMM J04433+0210 


4.5 


N5 


24.3 




MS 0440+02/0.19 





Note. — Tabic entries arc ordered by apparent 850-/j.m flux density. The most likely optical counterpart is listed first in column 3. The lensing 
amplifications were determined from the detailed mass models of the clusters. These were constructed using the LENSTOOL ray-tracing code 
on lensed features identified in high resolution optical images (Kncib et al. 1993). Uncertain redshifts are followed by a question mark. 



References. — a from Ivison et al. (1998) - b from Dressier et al. (1999) - c from Mellier et al. (1988) - d see also Gioia et al. (1998) - e sec 
also Dressier & Gunn (1992). 
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unobserved counterpart to SMMJ04433+0210). The red- 
shift distribution suggests that the majority of the sources 
are likely to lie at redshifts z = 1-3. The lens ampli- 
fication is not expected to significantly distort the N(z) 
distribution given the high redshifts of the bulk of the 
counterparts. In the bottom two panels of Fig. ^ we show 
for comparison the redshift distributions for the on-going 
Canada France Redshift Survey (CFRS, Lilly et al. 1999) 
and for the Hubble Deep Field (HDF, Hughes et al. 1998). 




1 2 5 



1 +z 

Fig. 5. — Cumulative redshift distribution (thick solid line), excluding 
cluster contamination, our two blank fields, and the unobserved source. 
Cumulative distributions for the three models discussed in Fig. 4 and in 
the text are also shown: Guiderdoni et al. (1998) Model E (dotted). Blain 
et al. (1999a) Gaussian model (dashed), and revised Gaussian model 
(solid; sec text). 



We next compare our observed redshift distribution with 
a few representative models formulated to fit the FIR back- 
ground and various infrared counts using simple analytic 
descriptions of the evolution of luminous FIR galaxies. 
We superimpose on our redshift distribution in Fig. ^ the 
Gaussian model (dashed line) from Blain et al. (1999a) 
and Model E (dotted line) from Guiderdoni et al. (1998). 
These models are composed of normal star forming galax- 
ies with an additional population of highly obscured galax- 
ies. 

In Fig. |5| we plot the cumulative distributions for our 
non-cluster data (also excluding the blank sources) and 
the models of Fig. ||[ When we perform the Kolmogorov- 
Smirnov (K-S) test on each model compared to the above 
data set, under the assumption that we have correctly 
identified the majority of the optical counterparts, we find 
low probabilities of 0.012 and 0.0004 for the Gaussian 
model and Model E, respectively; thus, both of these mod- 
els can be rejected for predicting a median redshift for the 
submm emission that is too high. However, when we do 
the same comparison with the reliable sample of four non- 
cluster counterparts (L1/L2, J1/J2, L3, and P4), we find 
higher probabilities of 0.53 and 0.11, respectively, neither 
of which are significantly inconsistent with the models. 

One way to reduce the mean redshift in the Blain et 
al. model while retaining the fits to the FIR background 
and counts is to change their dust emissivity spectral in- 
dex from 1.5 to 1.0. The 850-/im flux density of galaxies 
at z < 2 in the models is then increased as compared 
with galaxies at higher redshifts. Figures [l| and || show 
this modified model as a solid curve. The K-S test gives 
probabilities of 0.15 and 0.997 when this modified model is 
compared with the full and reliable samples, respectively. 






o ♦ 




b- — □ : 


- A 


- 




♦♦ * 




□ 




□ 


+ 






□ 


-i 





0.0 0.5 1.0 1.5 2.0 2.5 3.0 
Redshift 



Fig. 6. — Bolomctric luminosities derived from the lcnsing-correctcd 850- 
/^m fluxes assuming a dust temperature of 47 K (Arp 220), a dust emis- 
sivity spectral index of 1.0, and the lens amplifications listed in Table 2. 
Filled symbols correspond to secure optical identifications: the filled dia- 
monds are objects that show AGN activity, the filled star is our starburst 
pair, and the filled crosses arc the two cD galaxies. Open symbols cor- 
respond to uncertain optical identifications: the open cross is a cluster 
member, the open triangles arc 2 <C 1 field galaxies, and the open squares 
arc high redshift field galaxies. Note that the redshift range for J3 is in- 
dicated with a connecting line. 

Using the redshifts determined here, accurate cluster 
lens models, and a dust emissivity spectral index of 1.0, 
we can convert the apparent fluxes for the sources into in- 
trinsic bolometric luminosities assuming the Arp 220 dust 
temperature of 47 K (Klaas et al. 1997) and a qo = 0.5 cos- 
mology with Hq = 50 km s _1 Mpc~ . These luminosities 
are shown in Fig. || versus redshift and are characteristic 
of ultraluminous infrared galaxies (ULIGs). Since submil- 
limctcr observations of nearly all luminous infrared galax- 
ies have been reasonably fit by single temperature dust 
models with T = 30 - 50 K (Sanders & Mirabel 1996), 
the temperature dependence of the luminosity introduces 
a factor of only about 0.2-1.2 uncertainty in the above Lboi 
numbers. 

A very large fraction of local ULIGs show signatures of 
interactions and mergers. A relatively high fraction of dis- 
turbed or interacting counterparts were also uncovered in 
the optical identifications of the faint submm sample ana- 
lyzed here (Smail et al. 1998). Our spectroscopy has shown 
that although a small number of these interacting systems 
are simply projection effects (e.g. M1/M2), the majority 
are real; thus, we can infer that interactions remain an im- 
portant triggering mechanism for ultraluminous activity in 
the distant Universe. This conclusion supports the use of 
nearby ULIGs as templates to understand the evolution of 
these more distant systems. 

In Fig. ^ we show a more quantitative comparison of the 
properties of the local and distant ULIGs. We compare 
the observed optical to submm flux ratio versus redshift 
to that of the redshifted archetypical local ULIG Arp 220 
(solid curve). Interestingly, the four sources with firm 
identifications (L1/L2, J1/J2, Al835-cD, A2390-cD) all 
lie above the Arp 220 distribution in Fig. ^. However, a 
comparison of the properties of the reliable non-cluster 
counterparts (L1/L2, J1/J2, L3, and P4) indicates that 
there is considerable scatter in the Sj/Ssso ratio. For 
the non-cluster galaxies detected in the /-band, we find 
a median ratio of ~ 0.3 (S//Ss5o)Arp220 with a dispersion 
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of around a factor of four. The majority of our submm 
sources therefore appear to be emitting a slightly lower 
fraction of their luminosity in the optical relative to the 
submm as compared to that expected from an Arp 220- like 
source at their redshifts. We note that the blank sources 
would have a Sj/Sgso ratio of around 1 — 2 x 10~ 5 ; thus, 
if they are similar in their properties to our other sources, 
they would lie at redshifts z > 5. Alternatively, they could 
be more obscured systems at lower redshifts. 





A1835-cD 
O 


• 4<J 

▲ 3<j 




A2390— cD 
O 


o Foreground/Cluster. 




N1 \ 

o X 

HI 

o ° 






M6-M10 


\. L1-L2 

#M2 \, 






A K2 ^\ 







o 



^ 1 2 5 

1+z 

Fig. 7. — Ratio of /-band flux to 850-/im flux versus rcdshift. The line 
shows the relationship versus rcdshift expected for Arp 220. The galaxies 
show a large scatter in their relative Sj/Sgso ratios with the majority be- 
ing slightly fainter in the optical relative to their submm emission when 
compared to Arp 220 put at the same rcdshift. 

In terms of their bolometric luminosities, optical-to- 
FIR ratios, and morphologies, the galaxies selected in the 
submm in the distant Universe have very similar properties 
to local ULIGs. As regards the dominant energy source for 
the emission seen in the FIR background, i.e. AGN or star- 
burst, we find that at least three of the fourteen sources 



surveyed have counterparts with spectral features indica- 
tive of AGN activity. These objects may be part of the ob- 
scured AGN population predicted to be a major contribu- 
tor to the X-ray background at energies > 2 keV (Madau, 
Ghiscllini, & Fabian 1994; Comastri et al. 1995; Fabian 
et al. 1998; Gunn & Shanks 1999; Almaini, Lawrence & 
Boyle 1999). Deep radio, NIR, MIR, and high-resolution 
X-ray data of fields observed with SCUBA should provide 
us with improved position estimates and more information 
on the nature of the sources, making the task of following 
up the submm detections easier in the future. 

The determination of the redshift distribution of 
submm-selected galaxies is an important goal because it 
allows plausible models of the evolution of the volume 
emissivity of dust with redshift to be constrained more 
strongly than by the background and counts data alone. 
The first results on the redshift distribution presented here 
suggest that we can rule out models for galaxy and AGN 
evolution in which the bulk of the dust emission in the Uni- 
verse occurred at redshifts either below z ~ 1 or substan- 
tially higher than z ~ 3. This is in accord with conclusions 
reached by recent modelling of the FIR background (Dwek 
et al. 1998). Consequently, activity associated with these 
luminous dust-obscured sources occurs at a similar epoch 
to that seen for unobscured star-forming galaxies (Madau, 
Pozzetti, & Dickinson 1998; Steidel et al. 1999) and AGN 
(Boyle & Terlevich 1998). Our spectroscopic survey of the 
probable optical counterparts to submm sources is the first 
step towards the full reconstruction of the emission history 
of dusty galaxies in the Universe, whose contributions to 
star formation are now known to be at least as impor- 
tant at high redshift as the contributions from optically 
observed galaxies. 
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